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Entanglement Perturbation Theory for Infinite Quasi-ID Quantum Systems 
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We develop Entanglement Perturbation Theory (EPT) for infinite Quasi- ID quantum systems. 
The spin 1/2 Heisenberg chain with ferromagnetic nearest neighbor (NN) and antiferromagnetic 
next nearest neighbor (NNN) interactions with an easy-plane anisotropy is studied as a prototypical 
system. The obtained accurate phase diagram is compared with a recent prediction [Phys.Rev.B, 
81,094430(2010)] that dimer and Neel orders appear alternately as the XXZ anisotropy A approaches 
the isotropic limit A = 1. The first and second transitions (across dimer, Neel, and dimer phases) 
are detected with improved accuracy at A ~ 0.722 and 0.930. The third transition (from dimer 
to Neel phases), previously predicted to be at A « 0.98, is not detected at this A in our method, 
raising the possibility that the second Neel phase is absent. 

PACS numbers: 75.10.Pq , 75.10.Jm , 75.40.Mg 



In a recent article, Furukawa et al. [lj reported exotic 
Neel and dimer orders in a spin 1/2 Heisenberg chain 
described by the Hamiltonian, 
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characterized by a ferromagnetic NN exchange coupling 
Ji < and an antiferromagnetic NNN exchange coupling 
J 2 > 0, with an easy-plane anisotropy A < 1. Let us de- 
scribe the effect of the NNN exchange by /i = J2/J1. 
Figffl shows an equivalent coupled 2 spin-chains (zigzag 
chain) suited for EPT. The theoretical tools employed 
are the bosonization-Sine- Gordon (SG) theory [2j, the 
RG (renormalization group)-level spectroscopy analysis 
augmented by exact diagonalization 0, 0] , and numeri- 
cal RG analysis based on iTEBD (infinite time evolving 
block decimation) [5] , a variance of DMRG (density ma- 
trix renormalization group) method |6[ for infinite ID 
systems. The key finding is that the Neel order and 
dimer order alternate numerous times when approaching 
the ferromagnetic transition point (J2/I Ji|, A) = (1/4, 1) 
along the Lifshitz line which separates the character of 
the spin-spin correlation in the xy-plane from commen- 
surate to incommensurate. The reported Neel order is 
indeed counterintuitive since neither ferromagnetic NN 
nor antiferromagnetic NNN exchange couplings favor an- 
tiferromagnetism. 

When Js are both positive, the ground state phase di- 
agram of the model is well understood 0, 0, Q . Particu- 
larly, for small y, and with easy-plane anisotropy A < 1 , it 
is a gapless Tomonaga-Luttinger liquid (TLL). The gap- 
less TLL phase has instabilities toward ordered states 
Q. Neel order will be induced by easy-axis anisotropy 
A > 1, whereas dimer order appears when fj, is larger 
than some critical A-relatcd value. For instance, transi- 
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FIG. 1: (color online) A zigzag chain suitable for EPT. Cir- 
cles refer to spins. Long dashed lines (blue) indicate the 
NN coupling Ji, while short ones (orange) the NNN coupling 
■h = fJ-Ji ■ 



tion takes place when \i « 0.24 at A = 1. The quantum 
phase transition to these ordered phases can be under- 
stood within a bosonization-SG theory Q. 

Now in the case J\ < of our interest, the 
bosonization-SG theory of the TLL phase at J2 = is 
known exact, and the "infinite" alternation of energies 
of the Neel and dimer states was observed in the pa- 
rameter space < A < 1, and based on a perturbative 
treatment, from the Gaussian Hamiltonian, of the cosine 
term in the SG theory, the authors of [H Q claimed that 
the same should occur deep inside the frustrated regime 
~ 1/4, and thus the realization of the true alternating 
appearance of the Neel and dimer phases. This predic- 
tion was first checked by the level spectroscopy 0, 
showing 5 times of phase alternation between the two 
phases. This level spectroscopy analysis was, however, 
still based on the effective SG theory whose validity in the 
parameter regime |/x| ~ 1/4 is correct only in a pertur- 
bative sense such as neglecting a higher frequency term 
003(8-^7?). Thus, authors of [1| also carried out an in- 
dependent iTEBD calculation of the original model (1), 
confirming only one dimer to Neel transition at A ~ 0.72 



2 



and one Neel to dimer transition at A ~ 0.93. It was re- 
ported that iTEBD experienced a poor convergence for 
0.95 < A < 1 and 0.25 < \n\ < 0.35. 

The purpose of this letter is to study the controversial 
parameter regime in the above by EPT (entanglement 
perturbation theory), a recent many body method [Tol — 
1 1 31 ] . Our finding in this paper is summarized as follows: 
(1) EPT gives a converging, and hence exact phase dia- 
gram, including exact order parameters of Neel and dimer 
phases for the entire < A < 1 along the Lifshitz line, 
see Fig. 4 below. The dimer phase and Neel phase alter- 
nate only once, and for A > 0.93, only the dimer phase 
exists. This means two points. (2) First, the perturba- 
tive SG analysis when applied to the deeply frustrative 
regime \fi\ ~ 1/4, is missing some physics, and (3) sec- 
ond, iTEBD, even with a large entanglement x ~ 300, 
is not as accurate as EPT with \ ~ 80. This conclusion 
is not a surprise. Indeed, TEBD operates only on a few 
local spins, and therefore the deterioration of accuracy 
due to loss of information arising from successive Hilbert 
space truncation, typical of real-space RG, is expected 
worse than DMRG. 

The EPT algorithms EPT-gl and EPT-g2 both solve 
for ground state properties of systems with translational 
symmetry (a recent development of EPT handling inho- 
mogeneity will be reported elsewhere). EPT-gl is es- 
pecially suitable for infinite systems and is our starting 
point. It was previously applied to the Ji = case 
HI HI) and its superb accuracy is demonstrated in the 
first confirmation of a p rediction made by the confor- 
mal field theory [HI 17 1 on a long-range spin-spin cor- 



relation which can only be seen for thousands of lattice 
separation. It also confirmed yet another field theory 
prediction, by bosonaization, for the 4-spin correlation 
functions with better precision than DMRG fljj |. 

The central idea of EPT-gl is as follows. First is the 
system wave function in the matrix product state (MPS) 
lol 13- 23| representation 
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where is the local wave function matrix at site i, and 
N is the number of sites. cr;(=t, i) refers to local states. 
The dimension p of the square matrices , the entan- 
glement, controls the accuracy of the wave function \tp). 
In principle \ip) with optimized becomes exact with 
p — > oo. The matrices are optimized to maximize the 
variational energy: 
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value equation 
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where and Yi are symmetric matrices depending on 
{6}, and e is the eigenvalue. The eigenfunction {6} with 
the largest e describes the ground state wave function \tp). 



The density matrix K = e~@ H can also be written 
in matr ix p roduct form, the matrix product operator 
(MPO) |ldl. Il5j|. A generalization to quasi-lD naturally 
inherits this MPO of density matrix along with MPS for 
the wave function. But now, cf. Fig. 2, MPS is defined on 
a 2-spin composite (grey and black vertical two circles) 
and MPO is in a complex shell by shell (inner product di- 
rection) and layer by layer structure (vertical direction). 
We first deform a zigzag chain into two chains shown in 
FigfT] The Hamiltonian is rewritten as a summation of 
local bond Hamiltonians, 
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There are four types of bonds. Bond {1,2}, {3,4} etc 
form the first series. {2, 3}, {4, 5} etc belong to the 
second. These two series account for NN interactions. 
The third is {1, 3}, {2, 4}, {5, 7}, {6, 8} etc and the forth 
is {3, 5}, {4, 6}, {7, 9}, {8, 10} etc. They account for NNN 
interactions. The way to group is not unique, but it 
does not affect the result. It is immediately seen that 
[Hbondi, Hbondj] = if they are in the same series. This 
is the criteria to group the bonds. Small /3 safely sepa- 
rates series in the density matrix as follows 
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where /? — y [ll| . We then arrive at a generalized eigen- 



And for each series, we have 

(7) 

where repeating indexes follow Einstein summation con- 
vention. /, for instance for the XXX model, takes four 
2x2 matrices: 1, VJ3S X , y^S y , y^S z on a site. And g 
likewise on the other site of a bond. We use shell to de- 
note a series in (J7J. There are four shells in total. They 
are coupled in the direction of in and out of the paper in 
FiglU by the inner product of the operators on the same 
site. Also note that MPO contains two layers, composed 
of two vertically aligned spins, one in the top chain and 
the other in the bottom chain. We finally have an explicit 
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FIG. 2: Schematic figure of the MPO of the density matrix 
for a zigzag chain. Short dashed lines (orange) indicate NN 
interactions and long dashed lines (blue) NNN interactions. 
Solid lines refer to the inner prodoct. 




FIG. 3: Simplified schematic figure of the MPO of the density 
matrix for a zigzag chain. The two spins enclosed by an ellipse 
form an effective site. The dashed lines (orange) indicate the 
direct product and solid lines the inner product. 



expression of the whole density matrix as follows 



where 
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and r 2 likewise. a,/3 and 7 etc run from 1 to 4. They 
are the indexes involved in the direct product. If they 
are absorbed into an effective index, both T 1 and T 2 be- 
come 4-leg tensors. The two effective legs run from I to 
4 3 — 64. And the other two legs run from I to 4, respon- 
sible for the inner product with a new local wave function 
£ now defined on a 2-spin composite, denoted by an el- 
lipse in Fig 121 It is also noted that, reflecting the T 1 ' 2 
repetition structure, the local wave functions are also of 
two species £ 1,2 . The generalized eiegnvalue problem (4) 



is solved iteratively for the local wave functions, and af- 
ter convergence, energy etc can be calculated as before 

HE!. 



The calculation on zigzag chains is more expensive 
compared with a single spin chain because one now has 
to deal with the larger number of local states and larger 
indexes in MPO. Let us discuss the major time consump- 
tion. There are two. The first are generalized eigenvalue 
equations whose size is np 2 with n being the number of 
effective local states. The second is an eigenvalue decom- 
position of a transfer matrix in the horizontal direction 
in Fig. 3, a basis unit involving T 1 and T 2 and in total 
4 2-spin composites connected to r 1 ' 2 . Its size is mp 2 
where m is the size of the 3 bonds in the direct product 
between T 1 and T 2 . In contrast to n — 2 and m = 4 in 
a spin chain, they are 4 and 64 now. It is seen that the 
eigenvalue decomposition of mp 2 x mp 2 transfer matrix 
becomes especially large. 

Fortunately, there is a simple way to overcome this 
difficulty. The key point is in the original idea of EPT- 
gl, namely convert the Hamiltonian eigenvalue prob- 
lem to that for the density matrix K ee e~* H , noting 
e~P H — > 1 — (3H, leading to an MPO representation of 
the density matrix. A simple example is the product of 
two local density matrixes, (I — fiHij) ® (1 — PHkl) ~ 
I — fiHij — PHki + (3 2 Hij ® Hki , where the last term is 
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unimportant in the /3 — >■ limit. /3 = 10 _ 
used in this work like before [111 ]. Likewise, r*'? a /g'j> 
should contain such higher order terms in /3 which we 
can consistently discard. A rule to do so is in the fact 
that, in /„ and g a matrixes, a = 1 is identity matrix 
while other three terms are all order \J~f3. We should 
keep in mind that whenever there is a term of order 
there should be a companion order y/]3 term at the other 
end of the bond. This amounts to keep, among 64 index 
combinations in a/37, those terms with only one of a/37 
not being 1, namely, 111,112,121,211,113,131,311 and 
114, 141,411, there are total 10. 

To test the algorithm, we have calculated the dimeriza- 
tion d ee (S^i - Si+A ^ for fi = 0.5871 and A = 

where d is maximum when both Js are positive . The 
EPT result is 0.7906137 when p = 35 vs. 0.7906135 by 
DMRG |8j. By comparing with their results for /i = 2.5 
and A — 0, the most difficult case studied in [8] and prob- 
ably the same for EPT, it is found that EPT converges 
faster in entanglement. Since p-equivalent \ = 200 ~ 300 
is used in iTEBD the entanglement needed by EPT 
to achieve the comparable accuracy is much smaller than 
iTEBD. 

The main result of this letter is the phase diagram 
shown in Fig(4] where A approaches the isotropic point 
1 from below along the Lifshitz line. The Neel order 
parameter is {Sf) and the dimer order parameter D xy is 
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FIG. 4: The phase diagram along the line connecting 
(|M|,A) = (0.336,0.65) and (0.264,0.95), which approxi- 
mately traces the Lifshitz line. 



FIG. 5: Order parameters vs entanglement p at A = 0.722, 
one of the phase boundaries. Neel order wins the competition 
only at a large p. 



defined by 

D xy = (Sf.jiSf + Sf^Sf) — (SfS^ +1 + SfSf +1 ) (10) 

The main difference between our phase diagram and 
that of [1] is that our calculation converged solidly for 
all A up to 1, whereas iTEBD failed for 0.95 < A < 1. 
It also shows no more alternation of the Neel and dimer 
orders for A > 0.93, disproving the prediction of [H that 
there should be an infinite number of such transitions. 
The failed prediction, we argue, is most likely due to the 
failure of perturbative analysis of the Sine-Gordon theory 
deep in the frustrated regime ~ 1/4. Technically, the 
Lanczos diagonalization not only fails to converge due to 
highly degenerate nature for 0.995 < A < 1 as pointed 
out in 1], but also, we argue, yields not quite accurate 
energy levels for A > 0.95 so that the level spectroscopy 
analysis concludes incorrectly for this parameter regime. 
As for the nature of the Neel to dimer transition, our 
calculation, FigJH does not agree with [l] that it is first 
order. Rather it appears second order. 

While our calculation converged solidly, there is a clear 
evidence of a hard competition between the Neel and 
dimer orders near the transition point. In FigJSJ we have 
plotted the Neel and dimer order parameters as a func- 
tion of entanglement for A = 0.722. Dimer order dies out 
only at relatively large entanglement. The same happens 
at the other phase boundary A = 0.930 with even more 
subtle competition between the two orders since their 
magnitudes are both small. 

To summarize, we have developed EPT for quasi-lD 
strongly correlated quantum systems. There are two im- 
mediate applications of interest. First is the ID Hubbard 
model with long-range hopping, long range Coulomb in- 
teraction, orbital degeneracy and Hund-coupling in the 
pursuit of itinerant ferromagnetism and triplet supercon- 
ductivity. Second is a spin tube, namely the effective 



sites now contain the third spin, taking into account third 
nearest neighbor exchange interactions. Works are cur- 
rently under way, and will be reported elsewhere. 

Note added - After the completion of the paper, we 
have heard from the authors of [l| that the Lifshitz line 
along which we and they have done the calculations is 
actually not precise, and that their prediction of the 
Neel phase could still be observed close to the ferro- 
magnetic point, A = 1.0. We have explored A = 0.98 
and J2/IJ1I = 0.2552 (newly suggested point) to 0.2568 
(our and their previous calculation point). We observed 
a TLL (J2/IJ1I = 0.2552) to dimer order transition of 
the second order type which is indistinguishable from the 
BKT (Berezinskii-Kosterlitz-Thouless) type [H, Hj| due 
to an essential singularity [26j], the transition point being 
0.2560. We have also calculated the chiral order parame- 
ter S1XS2 in this parameter region to be strictly zero. In 
the meantime, we have observed a finite chiral order pa- 
rameter for the point A = 0.9 and J2/\Ji\ = 0.35 which 
is deeply in the chiral phase. We thus conclude that our 
main claim in the paper is unchanged. 

The numerical calculation utilized the RIKEN Inte- 
grated Cluster of Clusters at Advanced Center for Com- 
puting and Communication, RIKEN. We thank the au- 
thors of [l| for useful discussions. 
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